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Abstract. For random walks on networks (graphs), it is a theoretical challenge to 
explicitly determine the mean first-passage time (MFPT) between two nodes averaged 
over all pairs. In this paper, we study the MFPT of random walks in the famous T- 
graph, linking this important quantity to the resistance distance in electronic networks. 
We obtain an exact formula for the MFPT that is confirmed by extensive numerical 
calculations. The interesting quantity is derived through the recurrence relations 
resulting from the self-similar structure of the T-graph. The obtained closed-form 
expression exhibits that the MFPT approximately increases as a power-law function 
of the number of nodes, with the exponent lying between 1 and 2. Our research may 
shed light on the deeper understanding of random walks on the T-graph. 
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1. Introduction 

As the paradigmatic discrete-time realization of Brownian motion and diffusive 
processes, whose theory was formulated in the pioneering work by Einstein [1] and 
Smoluchowski [2] published in beginning of the last century, random walks have received 
a considerable amount of attention [HI HI E], uncovering their wide range of distinct 
applications. Thus far, random walks are still an active area of research [61 [TJ [HI El flO] - 
However, due to the complexity and variety of real media, theory of random walks 
on general graphs (networks) is not yet available. For this reason, studying random 
walks occurring on simple structures is a matter of exceptional importance. Fractal 
structures, in particular deterministic fractals, are valuable media in this content [7J[H], 
because their properties can be exactly studied. 

Among various deterministic fractals [HI [121 EH1 HU ESI [T6l [17] , the T- fractal (T- 
graph) [18] is a typical candidate for exactly solvable model, and a plethora of issues 
for random walks on this fractal have been studied [IHl EUJ EJJ E21 E31 EH EHJ EH EZ]- 
A simple analytical approach was proposed in [18], showing that random walks on the 
T-graph can be mapped onto diffusion on a comb structure. The asymptotic behavior 
of the moments of the first-passage time and survival probability for random walks 
on the T-graph was computed in [19]. Particularly, random walks performed on the 
T-graph with a single trap were extensively studied by many groups [191 EH EEl 127] . 
revealing some scaling relations and dominating behavior. Despite the fact that these 
investigations uncovered many unusual and exotic features of the T-graph, providing 
useful insight into understanding random walks on this fractal, they did not give a 
complete picture of random-walk dynamics on the T-fractal, since in these studies, only 
one special trap was considered. It was shown that in some real networks, any node may 
be looked as a trap, and the location of traps strongly affects the behavior of random 
walks [28]. Irrespective of its obvious importance and ubiquity, this issue has not been 
addressed for the T-graph. 

In this paper, using the connection between the random walks and electrical 
networks [29], we investigate random walks on the T-graph. The random walk process 
addressed here may be considered as the trapping problem with the perfect trap 
uniformly located at all nodes. We derive analytically an exact formula for the mean 
first-passage time (MFPT) averaged all pairs of nodes, which characterizes the efficiency 
of random walks on the T-graph. We show that the location of trap has no qualitative 
effect on the scaling for MFPT. We expect that our analytical method can be applied 
to some other deterministic media, and that our results can lead to deeper insights to 
random walks on the T-graph. 

2. Brief introduction to the T-graph 

The T-graph is built in an iterative way [181 [27]. We denote by T g (g > 0) the T-graph 
after g iterations. Note that henceforth we will also call the number of iterations as 
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Figure 1. (Color online) Iterative construction method of the T-graph. The next 
generation is obtained by performing the operation shown in the right-hand side of the 
arrow. 
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Figure 2. Illustration of the T-graph of generation 4. 



generation of the T-graph. The construction of the T-graph starts from (g = 0) two 
nodes connected by an edge, which corresponds to T . For g > 1, T g is obtained from 
T g -i by performing the operation illustrated in Fig. [IJ According to the construction 
algorithm, at each generation, the number of edges in the system increases by a factor 
of 3. Thus, we can easily obtain that at generation g, the total number of edges 
in T g is Eg = 3 9 . Since the T-graph is a tree, the total number of nodes in T g is 
N g = E g + 1 = 3 g + 1. Figure [2] shows schematically the structure of T 4 . 

The T-graph is a fractal with the fractal dimension and random walk dimension 
being df = p| ~ 1.585 and d w = |^| = 1 + df, respectively. The spectral dimension 
of the T-graph is d = 2df/d w — ~ 1.226. Note that for a general connected 
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graph, the spectral dimension governs the longtime behavior of a random walk on the 
graph at a deep level. For example, when a walker originating at a given node % of the 
graph, the probability Pu(t) for returning back to i, at long time t, obeys the relation 
Pu(t) ~ t~ d l 2 [30]. Again for instance, for an infinite graph with the spectral dimension 
d < 2, a walker starting from a given node will return back to the node almost surely 
over the course of time. The phenomenon of a random walk on such a graph is called 
"persistence". Since the spectral dimension d < 2 for the T-graph, a random walk on 
such a graph is persistent. 

To facilitate the description in what follows, we define the central node in Fig. [2] as 
the innermost node, and we call those nodes farthest from the central node as outermost 
nodes. Then, the T-graph can be alteratively constructed in another method, see Fig. [3j 
Given the generation g, T g+ i may be obtained by joining three copies of T g , denoted 
as T^, T g 2 \ and T g 3 \ respectively. In other words, to get T 9+1 one can merge together 
separate outermost nodes of the three replicas of T g . In the process of joining, the three 
outermost nodes are emerged into a single new node, which becomes the innermost node 
of T g+1 . 

3. Formulation of the problem 

In this section, we investigate a minimal model for discrete-time random walks of a 
particle on the T-graph T g . At each time step, the walker moves from its current location 
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to any of its nearest neighbors with equal probability. A key quantity characterizing 
such a random walk is the first-passage time (FPT), in terms of which many other 
quantities can be expressed. We are interested in the mean first-passage time (MFPT) 
between two distinct nodes, averaged over all pairs. 

It is well-known that the random walks addressed here can be described by the 
processes of Markov chains [31], the fundamental matrices corresponding to which can 
be used to express the FPT between any pair of nodes. However, the fundamental 
matrix method for computing the MFPT in T g requires computing the inversion of N g 
matrices of size (N g — 1) x (N g — 1), making it prohibitively difficult to calculate the 
quantity concerned for all but small networks. 

To bring down the high computational demands the fundamental matrix method 
makes, one can apply the method of the pseudoinverse of the Laplacian matrix [32] for 
the graph T g that random walks are performed on, which allows to compute the FPT 
between two arbitrary nodes directly from the network topology and requires inversion 
of a single N g x N g matrix. The elements Ifj of the Laplacian matrix L 9 for graph T g 
are defined as follows: Ifj = — 1 if the pair of two different nodes % and j is connected by 



a link, otherwise If, = 0; while l\f = k\ (i.e., degree of node %). Then, the pseudoinverse 



of the Laplacian matrix L„ is [33 
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where e g is the N g - dimensional "one" vector, i.e., e g = (1, 1, . . . , 1) T . 

Let Fij(g) denote the FPT for a walker in the T-graph T g , starting from node i, to 
arrive at node j for the first time. Since the graph T g is connected, the pseudoinverse 
matrix h g is well defined and the entries l\f can be applied to express the FPT Fij(g) 
as follows 



m=l 

where l 9 nm is the m entry of the diagonal of the Laplacian matrix ~L g as defined above. 
Then the sum, F tot (g), for the FPT between two nodes over all node pairs in graph T g 
reads 

F tot (g) = Y,J2 F M, (3) 

and the MFPT, (F) g , is 

The quantity of MFPT (F) g is very important since it measures the efficiency of 
the random walks on T g : The smaller the value (F) g , the higher the efficiency, and vice 
versa. Equations (T2J) and (jl]) show that the problem of calculating (F) g is reduced to 
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Figure 4. (Color online) Mean first-passage time (F) g as a function of the iteration 
g on a scmilogarithmic scale. The red squares are the numerical results obtained by 
direct calculation from Eqs. ([2|) and (j4|), while the full blue squares represent the exact 
values from Eq. (|23|) . both of which agree with each other. 



determining the elements of the pseudoinverse matrix ~L* , whose complexity is that of 
inverting an N g x N g matrix, it can be easily obtained by utilizing a standard software 
package, Mathematica 5.0. However, since N g increases exponentially with g, for large 
g, it becomes difficult to obtain (F) g through direct calculation using the pseudoinverse 
matrix, because of the limitations of time and computer memory. Therefore, one can 
compute directly the MFPT only for the first generations, see Fig. HI Fortunately, the 
particular construction of the T-graph and the connection [35j [36] between effective 
resistance and first-passage time allow us to calculate analytically the MFPT to obtain 
an explicit formula. Details will be given below. 

4. Rigorous solution to mean first-passage time 

In the following text, we will use the connection between the electrical networks 
and random walks to derive the closed-form expression for MFPT (F) g , avoiding the 
computation process for inverting any matrix. 

4-1. Relation for commute time and resistance distance between two nodes 

Given a graph, the underlying electrical network [29] can be obtained by replacing each 
edge of the original graph with a unit resistor. The effective resistance between any 
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node pair i and j is defined as the voltage when a unit current enters one node and 
leaves the other. It has been proved that the effective resistance is a distance measure, 
so the quantity is also called resistance distance [37] . Thus, we can view T g as a resistor 
network by considering all edges of T g to be unit resistors. Previous work [35], [36] showed 
that many problems on a resistor network are closed related to the classic random walks 
on the original graph. For example, the computation of effective resistance between two 
nodes in a resistor network can be expressed by the FPT between the two nodes on the 
corresponding network: For an arbitrary connected graph, the effective resistance Rij 
between a pair of nodes i and j is equal to (Fij + Fji) /(2E), where E is the total number 
of edges in the graph, and Fij is the expected time for a walker starting at node i to first 
reach j. In fact, F^ + Fji is the average time that a walker, originating at node i, will 
take to hit node j and go back to i, and it is often called commute time [HE] between i 
and j denoted as CV,-, i.e., CV,- = F^ + Fji. 

According to the close relation between commute time and effective resistance, we 
have that for T g , the effective resistance Rij(g) between a pair of nodes i and j is 
Rij(g) = Cij(g)/ (2 Eg), where Cij(g) = Fij(g) + Fji(g) is the commute time between i 
and j. Thus Cij(g) = Cji(g) = 2 E g Rij(g). Then, Eq. ([3]) can be rewritten as 

1 N g N g 

F **(9) = 2 E E C M = E * E E R ^3) , (5) 
and Eq. (J2J) can be recast as 

n^^EE^w^ («) 

Equation ([6]) shows that if we know how to determine the effective resistances, then we 
have the way to find the MFPT. 

For a general graph with order N, the complexity of resistance computation is that 
of inverting an N x N matrix, identical to that of calculating MFPT. However, for a 
treelike graph, the effective resistance between two nodes is exactly the usual shortest- 
path distance (also called geodesic distance) between the corresponding nodes. Since 
the T-graph has a treelike structure and it is self-similar at the same time, we can use 
these properties to find its geodesic distance, obtaining a rigorous expression. 

4-2. Exact expression for mean first-passage time 

After reducing the determination of MFPT to finding the average geodesic distance, the 
next step is to explicitly determine the latter quantity. To this end, we represent all 
the shortest-path distance of T g as a matrix [39], where the entry dij(g) is the shortest 
distance from node i to node j, which is the minimum length for the path connecting 
the two nodes. The maximum value D g of dij(g) is called the diameter of T g . Then, the 
average geodesic distance of T g is defined as the mean of shortest-path distances over 
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all pairs of nodes: 



S 

{L)d = N g (N g - l)/2 ' (7) 



where 



s g = Yl d M (8) 

denotes the sum of the geodesic distances between two nodes over all pairs. Note that 
in Eq. (jSJ), for two nodes i and j (i ^ j), we only count dij(g) or dji(g), not both. 

We continue by exhibiting the procedure of the determination of the total shortest- 
path distance and present the recurrence formula, which allows us to obtain S g+ i of 
the g +1 generation from S g of the g generation. According to the second construction 
method, we can easily see that the T-graph is self-similar. This obvious self-similar 
structure allows to calculate S g analytically. It is not difficult to find that the total 
shortest-path distance S g+ i satisfies the recursion relation 

S g+1 = 3S g + Q g , (9) 

where Q g is the sum over all shortest-path distances whose end-points are not in the 
same T g branch. The solution of Eq. (jHJ) is 

9-1 

S g = 3?S +Y l l&~ m ~ 1 nm]- ( 10 ) 

m=0 

Thus, all that is left to obtain S g is to compute fl m . 

The paths that contribute to Q g must all go through the innermost node (e.g., X 
in Fig. [3D of T 9+1 . To find fl g , we denote Qg'P as the sum of all shortest paths with 
end-points in and T^, respectively. Notice that Q g ,lB rules out the paths where 
end-point is the node shared by T g a ^ and T g . By symmetry, f2^' 2 = Q g )3 = fi^' 3 . Then, 
the total sum Q g is given by 

n g = &* + n 1 / + ni 3 = 3n 1 /. (n) 

In order to calculate the path length f^' 2 , we give the following notation. Let s g 
be the geodesic distance of a outermost node of T g jlQ], which is defined as the sum 
of geodesic distances between the outermost node and all nodes of T g , including the 
outermost node itself. Prior to determining s g , we first compute the diameter D g of 
T g . Actually, the diameter D g is equal to the path length between arbitrary pair of 
the outermost nodes previously belong to two different D g _i branches. Obviously, the 
following recursive relation holds: 

D g+1 = 2D g . (12) 

Considering the initial condition D = 1, Eq. (fl2j) is solved inductively to obtain 

D g = 2°. (13) 
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We now calculate the quantity s fl+1 . Let K denote an outermost node of T 3+1 , 
which is in the branch T^p (see Fig. [3]). By definition, s s+1 can be given by the sum 

s g+l = Y d Kj(9+^) 
jGT 9+ i 

= Y dK ^9 + 1) + Y dK ^9 + 1) + Y dK ^9 + 1) 

jeT« jeT ( 2 > J£T W 

= S , + 2 ^ d*;(«7+l), (14) 

where we have made use of X]j € T (2) j^x dK ^9 + -0 = SjgT (3) yj^ife + 1)' which is 
obvious from the symmetry. We denote the second term in Eq. (114jl by 2 ft, g . Then, 
= s g + 2h g . The quantity h g is evaluated as follows: 



h g = Y [ d Kx(g+ 1) + d X j(g+ 1) 



= ^ + (^-1)^, (15) 
where dxx(g + 1) = was used. Hence, Eq. ([Hj) becomes 

Sff+1 = 3 Sg + 2(iV s -l)D ff . (16) 
Using iV 5 = 3 3 + 1, -D g = 2 s , and s = 1, Eq. fflBT) is resolved by induction 

s 9 = 3 9 - 1 (2 3+1 + l). (17) 

With above obtained results, we can determine ^ 9 ' 2 , which can be expressed in 
terms of the previously explicitly determined quantities. By definition, Q g ' 2 is given by 
the sum 



n J' 2 = Y d ^9 + 1) 



Y [d iX (g + l) + d Xj (g+l) 



j&ff\j^X 

= (Ng ~l)Y ^9 + 1) + (Ng ~l)Y ^9 + 1) 
= 2{Ng-l)Sg 

= 2 x 3 29 ~ 1 (2 S+1 + 1), (18) 

where we have used the equivalence relation ^ ieT (i) { , x d iX (g+l) = X^ e T (2) j^x d xj(9+ 
!)• 

Inserting Eq. (TT81) into Eq. ( TTTl) . we have 

0, = 2 x 3 29 (2 9+1 + 1) . (19) 
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(L) 9 = ,Z,On,J ■ (21) 



F tot (s) =2£ 9 5 3 = -x 9 9 (4 x 6 9 + 5 x 3 9 + 6) (22) 



Substituting Eq. (TT9~j) into Eq. (fill and using the initial value So = 1, we can obtain 
the exact expression for the total distance 

5 s = ^xr 1 (4x6 s + 5 x 3 9 + 6). (20) 

Then the analytic expression for mean geodesic distance can be obtained as 

8 x 6 s + 10 x 3 9 + 12 
15 x (39 + 1) 

Since for T ff , the shortest-path distance is equivalent to the resistance distance, i.e., 
dij(g) = Rij(g), according to the above obtained results shown in Eqs. (J20l) and (l2Pj) . 
we can easily get that the total of FPT between all N g (N g — 1) pairs of nodes and the 
MFPT are 

2 
15 
and 

2 9 3 9 

(F)„ = — 9 - = EJL) V = — -(8 x 6 9 + 10 x 3 9 + 12) , (23) 

\ ig Ng g\ 19 15(39 + >> v > 

respectively. 

We have checked our analytical formula against numerical values, see Fig. |H For 
the range of 1 < g < 8, the values obtained from Eq. fl23|) are in complete agreement 
with those numerical results on the basis of the direct calculation through the method 
of pseudoinverse matrix discussed in section [3j This agreement serves as an independent 
test of our theoretical formula. 

Note that the result provided in Eq. (1231 is consistent with the recent proposal by 
Condamin et al. [JT] for a very general scaling form for FPT, (F), of a random walk as 
a function of the distance, L, between the origin and the trap location. Condamin et al. 
proved that, for complex invariant networks with d w > df, the leading behavior of FPT 
behaves as (F) ~ NL dw ~ df . For the T-graph, d w — df = 1, according to the conclusion 
of Condamin et al, one has (F) g ~ N g (L) g as found in Eq. ( 1231) . 

Next we will express the MFPT (F) g as a function of network order N g , in order 
to obtain the scaling between these two quantities. From the relation N g — 3 s + 1, we 
have 3 9 = N g — 1 and g = log 3 (N g — 1). Thereby, we can rewrite Eq. (|23l as 

(F) g = 2( ^~ g 1} (4(iV g - l) 1+log 3 2 + 5(N g - 1) + 6) , (24) 

which provides explicit dependence of {F) g on network order N g . In the infinite network 
order, namely N g — >• oo, 

(F) g ^^(N g y^ = ^(N g )^. (25) 

Thus, for the large T-graph, the MFPT grows as a power law function of network order 
with the exponent larger than 1 and less than 2, implying that the MFPT increases 
superlinearly with the number of network nodes. 

Equation (1251) encodes the speed of the random walks on T g , which may be 
quantified by the coverage C g (t), standing for the mean number of different nodes 
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visited by a walker at time t, averaged for distinct walks initially starting from different 
sources [121 E2]. As t — > oo, the leading asymptotic behavior for C g (t) is 

C g {t)~t*'\ (26) 

which may be established from the following heuristic argument [HI B5] . Since the 
random walker will visit all nodes of T g with probability one in the longtime limit, 
irrespective of its starting location, N g plays the role of C g (t), and (F) g is like t (step 
number of the walker). Then, inverting the relation in Eq. (126]) leads to asymptotic 
scaling (F) g ~ (N g ) 2 ^ d , as found above and given by Eq. ( |25|) . which also provides the 
exact coefficient of proportionality. 

Notice that in Ref. [27] Agliari studied the trapping problem with an immobile 
trap located at the innermost (central) node, and showed that in the asymptotic limit, 
the average trapping time (ATT) to first reach the trap, averaged over all nodes except 
the absorbing node itself, exhibiting a similar behavior (the same exponent) as that of 
Eq. (f25|) . but has a different prefactor. Note that, due to the symmetry, the ATT in 
this case corresponds to the mean time for random walks with a fixed trap situated at 
an outmost node in the T-graph of generation g — 1 [27] . That is to say, random walks 
of these two extreme cases show almost the same behavior of ATT. On the other hand, 
the random walks discussed here may be considered as the trapping issue with the trap 
uniformly distributed throughout all nodes on the T-graph. Thus, we can conclude that 
the trap's location has no qualitative impact on the leading behavior of the ATT for 
random walks on the T-graph, which is in sharp contrast to trapping problem on some 
scale- free networks [28], where the ATT depends on the trap position: it is much smaller 
for the case when trap is located at a large-degree node than that for the case when the 
trap is uniformly placed. 

5. Conclusions 

In this paper, we have studied the standard discrete-time random walks on the T-graph. 
By making use of the link between the electric networks and random walks, we have 
studied analytically the MFPT averaged over all pairs of nodes, based on the recursive 
relations derived from the self-similar structure of the T-graph. We have determined 
the rigorous solution for the MFPT, which in the large limit of graph order, increases 
algebraically with the number of nodes. We showed that trap location has little influence 
on the scaling of MFPT for random walks on the T-graph. We expect that by providing 
a paradigm for computing the MFPT, our analytical technique could guide and shed 
light on related studies of random walks on other deterministic media, especially to some 
scale- free fractals jl6j Wf\ that has received much recent contention [%8l 1^9] . We also 
expect that our work, especially the exact solution, could prompt the studies of diffusion 
process on random fractals, other stochastic graphs, and even complex networks [50J, 
by giving a guide to and a test of approximate methods for random media. 
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